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f-*^ ' Abstract. We investigate stresses and particle motion during the start up of 

flow in a colloidal dispersion close to arrest into a glassy state. A combination of 
^ , molecular dynamics simulation, mode coupling theory and confocal microscopy 

el experiment is used to investigate the origins of the widely observed stress 

' overshoot and (previously not reported) super-diffusive motion in the transient 

dynamics. A link between the macro-rheological stress versus strain curves and the 
microscopic particle motion is established. Negative correlations in the transient 
, auto-correlation function of the potential stresses are found responsible for both 

O . phenomena, and arise even for homogeneous flows and almost Gaussian particle 

vJ 1 displacements. 
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lO 1- Introduction 

(N 

0^ ' Relaxation processes in highly viscous liquids are strongly affected by the application 

^ ] , of an external shear field. Even at small shear rates, the shear viscosity can decrease 

t^^ ' by orders of magnitude while the structure exhibits only small changes compared to 

^^ . equilibrium. This phenomenon of shear thinning has been observed in experiments 

of various systems, such as colloids and polymers, as well as in computer simulations 
of model systems [l| . The transport coefficients measured in quiescent equilibrium or 
linear response markedly differ from those in steady-shear states. This raises the 
question how the system's microscopic dynamics evolves in response to a sudden 
^ ' change in the externally applied shear. This concerns the transient dynamics, as, for 

example, displayed in time-dependent transport coefficients or waiting-time dependent 
two-time correlation functions. Here, we are concerned with a suddenly commencing 
shear flow, imposed upon a quiescent, equilibrated system. While the measurement 
of macroscopic quantities like the shear stress a is standard rheology procedure, the 
detailed study of microscopic dynamics has been a recent contribution to the field 
from both computer simulation [2, y, \4, la, lo, LZl, M, M ^^id direct-imaging techniques 
for colloidal suspensions, namely confocal microscopy [lO|. 

It is well known that after starting a steady shear fiow, stresses develop in the 
liquid. At short times after switching on (corresponding to small overall strains 
7 — jt), they increase almost linearly and at long times saturate at a constant level, the 
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stress corresponding to the sheared steady state. In many cases, the increase in stress 
first continues beyond this steady-state value before relaxing back at larger times. 
A nonmonotonic cr(i)-versus-i curve results, with a maximum at intermediate strains, 



the 'stress overshoot'. Such overshoots have been seen in polymeric liquids [11|.|12|.|13||. 
metallic glasses [IJ, ll5| , and in simulations [6|, |l6|, [iTI, UM, UM ■ Despite its abundance 
among different systems, no generally accepted physical picture has emerged. Linking 
this macroscopically observed feature to the underlying microscopic dynamics can be 
achieved by studying the transient (microscopic) correlation functions. We pursue 
this by combining confocal microscopy, molecular-dynamics simulation, and a mode- 
coupling-theory approach. Surprisingly, we find that a previously not reported, super- 
diffusive particle motion is characteristic for the dynamics at strain values where the 
stress overshoot appears. We connect both phenomena to negative portions in the 
transient auto-correlation function of the microscopic stresses. 



The mode-coupling theory (MCT) for colloidal rheology [20|, l2l|, 123, l23|, l2J is 
a recent microscopic framework from which constitutive equations can be derived 
that capture the interplay of slow dynamics and shear thinning. The theory can 
thus address the above question on the relation between microscopic dynamics and 
macroscopic rheology. Its underlying physical picture is a combination of the cage 
effect, leading to slow dynamics, and of the advection of long- wavelength fluctuations 
to short wavelength, leading to a shear-induced breaking of cages 25|, |26|, |27l |28[ . The 



MCT equations can be formulated from an integration-through-transient formalism 



developed by Fuchs and Gates [20|, |2l|. In this formalism, transient correlation 
functions of a special kind (equilibrium averages over fluctuations that evolve according 
to the nonequilibrium dynamics) play a central role. The direct measurement of such 
correlations, which is possible under shear startup, provides a direct link between 
theory, simulation and experiment. This allows for a detailed test of some of the 
assumptions implicit in the theory. 

The paper is structured as follows: we will in section [2] discuss the different 
techniques used in this study, viz. the simulation model, experimental setup, and 
theoretical framework. Section [3] presents and compares the main results for 
both stress buildup and waiting-time dependent mean-squared displacements after 
switching on shear flow. Finally, section |4] concludes. 

2. Methods and Techniques 

2.1. Simulation model and details of the simulation 

Molecular dynamics computer simulations have been done for a simple model of a 
binary AB mixture of charged colloids. The interactions between the particles are 
modeled by a Yukawa potential, 

exp(-Ka^(r-rfa/3)) fl A T3 n\ 

Uap ^ eapdap a,fi = A,B, (1) 

which is truncated at a cut-off distance r"'^, defined by Uap{r"'^) = lO^^eAA- The 
"particle diameters" are set to d = c^aa = 1-0, ^bb = 1-2 d and ^ab = 1-1 d, the energy 
parameters to e = eAA = 1-0, ebb = 2.0 e, eab = 1.4 e, and the screening parameters 
to kaa = '*BB = '^AB = &/d. The choice of these parameters ensures that, at the 
density g = 0.675 mj^/d\i^^ considered in this work, no problems with crystallization 
or phase-separation occur, at least in the temperature range under consideration. The 
masses of the particles are set to unity, i.e. m = toa = m-Q = 1.0. 
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The simulations were done for a 50:50 mixture of A^ = 2Np^ — 2Nb = 1600 
particles, placed in a cubic simulation box of linear size L = 13.3 d. For the sheared 
system, we chose the x direction as the direction of shear and the y and z direction as 
the gradient and vorticity direction, respectively. Shear was imposed onto the system 
via modified peri odic boundary conditions, the so-called Lees-Edwards boundary 
conditions 29|,l30|. Here, a particle that moves out of the simulation box in y direction 



is subject to a displacement in x direction due to constant velocities Us^x and —Us^x of 
the image cells above and below the simulation cell, respectively. For the glassforming 
Yukawa system considered in this work, the application of Lees-Edwards boundary 
conditions leads to a linear shear profile in the steady state regime, fs,£c(2/) = 'iiy~L/2) 
with the shear rate 7 = Us^x/L. In recent simulation studies, Lees-Edwards boundary 
conditions have been used in conjunction with the so-called Sllod equations of motion 
[30| to enforce the formation of a linear shear profile. For our purpose, the use of the 
Sllod equations is not appropriate since we are mainly interested in the study of the 
transient dynamics, i.e. in the time regime before the steady state regime is reached. 
In this case, one may expect that the emergence of non-linear shear profiles strongly 
affects the response of the system to the external shear field. Thus, the Sllod equations 
would modify the transient dynamics in an artificial manner. However, as we shall 
see below, an almost linear profile is built up long before the steady state regime is 
approached. So our simulations indicate a posteriori that it does not matter whether 
or not the Sllod equations are used. We also note that non-linear shear profiles may 
occur as a consequence of a hydrodynamic instability, leading to the formation of shear 
bands. However, in the glassforming Yukawa system of this work, such phenomena 
are not observed, at least for the considered range of shear rates. 

Both for the equilibrium simulations and the simulations under shear, the system 
was coupled to a dissipative particle dynamics (DPD) thermostat (see Ref. [31| and 
references therein) . The DPD equations of motion are given by 

^^ = ^' P^=E[F..+FS+FS] (2) 

with Ti and p; the position and momenta of a particle i [i ^ l,...,iV). In Eq. ^ 
Fij = —Vuij denotes the conservative force between a particle i and a particle j due 
to the interaction potential defined by Eq. ([T]). To provide the thermostatting of the 
system a dissipative force Ff^ and a random force F^ are added in Eq. ^. 
The dissipative force is defined by [3^ 

Ff^=-Cw^r,,) (f,j.v„)f,, (3) 

with C being a friction coefficient, v^j = v,; — Wj the relative velocity between particle 
i and j, Vij the unit vector of the vector r,;j = r,; — Vj, and rij = jr^ — r^l the distance 
between particle i and j. The function w{rij) is defined by w{rij) = ^1 — r^j /re 
for rij < r^^^ = 1.25 d and w{rij) = otherwise. Thus, the force F^ describes a 
frictional force due to the interaction between neighboring particle pairs. The use 
of the relative velocities between neighboring particles in Q is crucial to obtain the 
correct behaviour on hydrodynamic scales. It ensures that the DPD thermostat is 
Galilean invariant and that momentum is locally conserved. For the friction coefficient 
C we chose the value C = 12. With this value, the microscopic properties are close 
to those of a purely Newtonian dynamics with C = 0. The microscopic dynamics 
can be significantly changed by using high values of C, approaching the limit of an 
overdamped stochastic dynamics where inertia effects can be effectively neglected. In 
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a forthcoming publication [33| . we present simulations with ^ = 1200 and compare 
them to those with C = 12 that are shown in the present paper. These simulations 
indicate that qualitatively all the essential features of the transient dynamics do not 
depend on the choice of (. 

The random force in ([2]) is given by F^- ~ ^/2k^rC,w{rij) 9ij Vij where 9ij = 9ji 
are uniform random numbers with zero mean and unit variance. The amplitude of 
the random force, ^/2kBT(, is chosen in accordance with the fluctuation-dissipation 
theorem. 

The equations of motion were integrated by a generalized form of the velocity 



Verlet algorithm that has been recently proposed by Peters [3J| . For the time step of 
the integration we used St = 0.0083 r (with the time unit t = ^JmSFje). First, the 
samples were fully equilibrated in the temperature range 1.0 > T > 0.14, performing 
at least 30 independent runs at each temperature. Starting from fully equilibrated 
samples, production runs at various shear rates were performed. These runs were 
sufficiently long to reach the steady state regime. 

At the lowest temperature T = 0.14, 250 independent runs were done, each of 
them over 40 million time steps. This relatively large effort was necessary for an 
accurate determination of the shear stress, defined by 



^i: 



L3 



j>i 



(4) 



The shear stress (axy) is a collective quantity and thus lacks the self-averaging property 
of one-particle quantities. Therefore, a relatively large number of independent runs 
have to be made to obtain (cr^y) with a reasonable accuracy. Below we also present 
simulation results for the mean squared displacement (MSD) of a tagged particle, 

(r2(t)) = ((rtag,a(t) -nag,Q(0))^), a^x,y,z, (5) 

where rtag^ait) is the a component of the position of the tagged particle at time t (an 
average is performed over all particles of the same species). The MSD is an example of 
a one-particle quantity. From the MSD. the self-diffusion constant D can be calculated 
via the Einstein relation, 

i^.limilM. (6) 

t^oo 2t 

Note that for the sheared system, the MSD as well as the self-diffusion constants are 
anisotropic and depend on the considered Cartesian direction (see below). 

Another quantity, that we use to characterize the difference between sheared and 
unsheared systems, is the so-called non-Gaussian parameter [35l. l36l|. 

a2 (t) is the coefficient of first correction term to the Gaussian approximation of the 
incoherent intermediate scattering function. Note that 0:2 (i) vanishes in the diffusive 
long-time limit when the Einstein relation ^ holds. 

2.2. Experimental techniques 

2.2.1. Samples We used polymethylmethacrylate (PMMA) colloids fluorescently 
labeled with nitrobenzoxadiazole (NBD) and dispersed in a mixture of cycloheptyl 



From Equilibrium to Steady State 5 

bromide and cis-dccalin that closely matches the density and refractive index of the 
colloids. Since the colloids acquire a small charge in this solvent mixture, we added 



4 mM of tetrabutylammonumchloride in order to screen the charges [37[ . Such a system 
shows nearly hard-sphere (HS) behaviour, where the volume fraction cj) = {7r/6)n(P is 
the only thermodynamic control parameter. 

The particle diameter d = 1680 ± 4nm was determined by static light scattering 
on a very dilute colloidal suspension (0 ~ 10^'*). The diameter was also independently 
estimated from the position of the first peak of the radial distribution function obtained 
by confocal microscopy [38|, l39[, yielding d = 1690 nm. With the very dilute sample 



we also performed dynamic light scattering. This confirmed the determined particle 
diameter and, in addition, allowed us to deduce the relative polydispersity in size, 
0.062, from the angular dependence of the diffusion coefficient [40|. 

The colloid volume fraction of the stock solution was calibrated by drying. A drop 
of the suspension was weighed and allowed to dry in a vacuum oven until all the solvent 
was evaporated. The weight fraction was then calculated as (j)^, = '«dry/'7itotai, where 
mdry is the mass of the dried sample and TTitotai the original mass before drying. The 
colloid volume fraction was estimated from the weight fraction as = acj)^ , where the 
factor a takes into account the contribution of the PHSA hairs, which are collapsed 



in the dried state. Following [4lJ we assume a = 1.04. Our sample was prepared at 
(j) = 0.57. The volume fraction of the sample was also determined by direct imaging 
and found to be consistent. 

2.2.2. Shear Cell Shear was applied to the sample by means of a shear cell designed 
and built in our lab. It represents an improved version of cells previously used in 
light scattering and microscopy experiments 42, l43l . l44l | . For the present microscopy 



experiments, a specially designed microscope stage is used to mount the shear cell in 
the microscope. The cell consists of two parallel glass coverslips (thickness 170 /im) 
whose separation can be continuously varied within the range 300 ... 1000 /im by 
vertically shifting the upper plate, which is sitting on three metal spheres. Adjusting 
the position of the metal spheres also ensures that the two plates are parallel; this 



is verified as described before [45|. In the present experiment we used a gap width 
(320 ^m) close to the minimum value in order to achieve maximum strain. The bottom 
plate of the cell is driven by a piezoelectric actuator (PI Instruments, model P-841.60) 
which can provide a maximum stroke of 90 yum. The other end of the bottom plate 
is connected to a pivoted lever driving the upper plate in the direction opposite to 



the bottom plate [421 . The vertical position of the pivot moves the vertical position 
of the stationary plane inside the sample and changes the maximum displacement of 
the top plate. In this experiment the maximum displacement of the top plate was 
630 ^m, as measured by a displacement sensor (Lion Precision, model ECL100-U8A, 
displacement range 2 mm) which was calibrated prior to the measurement using the 
piezoelectric actuator displacement, which is known with nanometer accuracy. With 
this arrangement the maximum obtainable strain is 7 = 2.25. 

A constant shear rate is achieved by driving the piezoactuator with a linear voltage 
ramp which covers the whole displacement range. The slope of the ramp determines 
the shear rate 7. The shear rate is also experimentally measured from the velocity 
profile across the gap, as explained below. 

In order to minimize solvent evaporation during the experiments, the two glass 
plates are surrounded by teflon sheets which touch each other and thus seal the sample 
environment. Due to the very low friction of teflon, no change in the plate displacement 
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Figure 1. Confocal microscopy image of the sample at switch on of shear. 
The axes x and z indicate respectively the velocity and vorticity direction. The 
gradient direction is perpendicular to the plane of the image. 



due to the contact of the teflon sheets is observed, as verified during cahbration. 

In recent investigations of the flow of glassy coUoidal dispersions wafl shp was 
observed when smooth surfaces were used to apply shear to the samples. It was shown 
[10| that slip can be prevented when the surfaces are coated with colloidal particles 
similar to those used for sample preparation. Although our volume fraction is below 
the glass transition, we coated the coverslips with a layer of polydisperse colloids. The 
velocity proflle across the gap indeed shows no indication of wall slip (see below). 



2.2.3. Confocal Microscopy and Image Analysis Confocal Microscopy experiments 
were performed with a fast-scanning VT-Eye confocal microscope (Visitech 
International) mounted on a Nikon TE2000-U inverted microscope. A Nikon Plan Apo 
VC 100 X oil immersion objective was used for all measurements. Two-dimensional 
images of the samples were recorded at a depth of 15 //m inside the sample (except 
for the determination of the velocity profile, see below), in order to avoid boundary 
effects and to retain good quality images. Images have 512 x 512 pixels, corresponding 
to an area of 57 x 57 /im^. 

In a typical experiment 500 images were recorded starting simultaneously with 
the application of shear and ending when the maximum strain was reached. Image 
acquisition and application of shear were synchronized by starting the voltage ramp, 
which controls the displacement of the plates, with a trigger provided by the confocal 
microscope. The response of the piezoactuator to the voltage (/iS delay) can be 
considered instantaneous on the time scale of our experiments. 

Figure [1] shows a typical confocal image acquired in experiments. The axes x 
and z indicate the velocity and vorticity direction. The gradient direction, y, is 
perpendicular to the plane of the image. Particle coordinates and particle trajectories 
were extracted from images using standard routines [46|. The shear rate applied here 
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Figure 2. Histogram of particle displacements between consecutive images along 
the velocity x (o) and vorticity z (I) direction. 



(7 = 0.045 s~^) is sufficiently small so that particles do not move very far between 
two consecutive frames even in the x direction. It is thus not necessary to remove 
the affine motions before tracking [47| . This is supported by the histograms of the 



frame by frame displacements in x and z directions as determined during a step rate 
experiment (Fig. [5]); both histograms lie well within the displacement axis, which 
covers the maximum detectable range (about ±1 fim). 

Particle trajectories were used to calculate mean-squared displacements (MSD) 
and non-Gaussian parameters (0^2) in the vorticity direction z, according to equations 
[5]and[7]respectively. Only trajectories starting at commencement of shear are included 
in the MSD and a2 calculations in order not to loose information on the time 
dependence of the MSD and a2 after switching on the shear field. This implies that 
in each experiment only a limited number of particles are contributing to the average 
value of the two quantities; initially of the order of 1000 and the number decreasing 
with time as particles leave the field of view due to motion in or out of the observation 
plane. In order to improve statistics, MSDs and 012 extracted from several, typically 
20, experiments were averaged after checking reproducibility. 

The velocity profile across the gap was determined immediately after the 
application of shear by recording image series at different depths into the sample. The 
depth was increased in steps of 10 /im, from the first layer of particles above the coating 
up to 100 fim into the sample. Deeper in the sample image quality did not allow for a 
quantitative analysis. Displacements between consecutive images were determined by 
maximizing the correlation between consecutive images j47| . The obtained velocity 
profile (Fig. [9]) shows a linear dependence with a shear rate 7 = 0.0458 s~^, in good 
agreement with the nominal value (for a detailed discussion sec below). 



2. 3. Mode coupling theory for glassforming liquids under shear 



Recent advances [20|, |21|, |23|, |2J] have generalized the mode coupling theory of the glass 
transition to describe dense colloidal suspensions in a fiowing solvent, characterized 
by a non- vanishing velocity gradient tensor of arbitrary time dependence K(t). The 
resulting non-linear theory provides a route to calculating non-equilibrium averages 
of functions of the particle coordinates. In particular, a closed microscopic expression 
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for the shear-stress in terms of the fuU flow history of the system, i.e., a eonstitutive 
equation relating stress and strain, may be obtained. The starting point of the theory 
is the many-body Smoluchowski equation [48j 

dt^it) = n{t)^it) = ^ a, • [Doid^ - PFi) - K{t) ■ r,;] *(t) , (8) 

i 

where Da is the bare diffusion coefficient and /3 = l/ksT. Equation jS]) describes 
the evolution of the probability distribution function of particle positions ^(i) in the 
overdampcd limit. This formulation of the theory omits hydrodynamic interactions 
and assumes a prescribed, spatially constant /«(i), thus excluding confinement 
effects and inhomogeneous flows such as shear-banded and shear-localized states. 
Although the assumption of translationally invariant flow is clearly an idealization, 
the simulation and experimental results presented in this work suggest that Eq.® is 
a reasonable starting point, at least for dense fluid states under shear flow. In the case 
of a constant shear flow 7 along the a;-direction and with gradient in the j/-direction 
the velocity gradient tensor is given by Kij ~ j SixSr 



The integration through transients formalism [20|, |2l|, |23|, |24| provides a formal 
solution to Eq.® by integrating over the entire flow history. In the case of shear, the 
following formal result is obtained for the time-dependent distribution function 

*(t) = *e + / dt' Jit') vl/ea,,e{'' ''■''' ^^^ , (9) 



where a^y = — ^^ Fi^riy is the potential part of the microscopic stress tensor and 
fit is the adjoint Smoluchowski operator. The time-ordered exponential function 
e_ arises because fl^t) does not commute with itself for different times. The 
assumption of an equilibrium Boltzmann distribution ^e in the inflnite past allows 
general nonequilibrium averages to be expressed in terms of averages taken with the 
equilibrium distribution function. Equation ^ is an operator expression, to be used 
with the understanding that quantities to be averaged are placed to the right before 
integration over particle coordinates. Within linear response the macroscopic shear 
stress a{t) is given in terms of the shear modulus by the familiar Green-Kubo formula 
of equilibrium statistical mechanics. Using ([9]) to calculate the average of the potential 
part of the macroscopic stress, viz. (T^i-y/y, where V is the system volume, yields a 



non- linear generalized Green-Kubo relation [20|, |2l|, [23|, [24 1 



*^ J —CO 

where (•) represents an average over the equilibrium distribution. Note that the 
flow history 7(i') appears both explicitly in Eq. p^ (recovering linear response) 
and implicitly (nonlinearly) through the time- evolution of the stress-stress correlation 
function. In the present situation of start-up shear, 7(t') vanishes for all t' < 0, so 
that the integration is only performed for t' > 0, where the Smoluchowski operator 
has no explicit time dependence. Defining the generalized dynamical shear modulus 
by G{t) = {axy{t)(Txy{Q)) /V as a non-linear function of 7, we have 

<Jxy{t)^^ fdt'G{t'). (11) 

Jo 
Within the mode-coupling approach G{t) is approximated by projecting the dynamics 
onto density-pair modes corresponding to all possible wave- vector pairs and directions. 
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The considerable numerical complexity of the resulting equation can be significantly 
reduced by employing the following isotropic approximation for the modulus [22| . 

where Sk is the static structure factor. The affine solvent flow enters via the 
time-dependent wavevector k(t) describing the advection of density fluctuations for 
wave number k to smaller wavelengths. The exact anisotropic advection k(t) = 
■\/k^ + 2kxky'yt + k'^jH'^ needs to be approximated by an isotropic wavevector in 
the mode-coupling fc-integrals, k{t) « fc-^/l + (7t/7c)^/3, in order to reduce the 
computational resources required. Pre-averaging over spatial directions introduces 
an additional source of error when compared to full solution of the anisotropic 
MCT equations. To first order, we expect, this discrepancy can be compensated by 
introducing a characteristic strain parameter, 7c, to which we assign the value jc = 0.1. 
Evaluation of (|12p requires the transient intermediate density correlator, defined by 
^k{t) = {g'^cxp['[lH]gk(^_t)) /NSk- The appearance of an advected wavevector in 
this definition removes the trivial decorrelation of density fiuctuations arising from 
purely affine flow. As a consequence of the equilibrium averaging, ^k (t) only contains 
information regarding the strain accumulated between t = and later time t. Mori- 
Zwanzig type projection operator manipulations yield an equation of motion for the 
correlator [231 



-<f>(t) + r, (^$^(t) + ^ dt' m,{t - t')-^<i>,{t')^ = 0, (13) 



dt 

where Tq ~ q^D^j Sq is the initial decay rate. Mode-coupling approximations provide 
an explicit form for the memory function 

mq{t) = dk^^--^[q-kck(t) +q-pCp(t)] [q • kc^ + q • pcp], (14) 

where n is the number density and Cq is the equilibrium direct correlation function. 
ncq = 1 — l/5<j. It should be noted that the colloid-colloid interaction potential and 
thermodynamic statepoint enter purely via the equilbrium structure factor. In the 
absence of flow the two vertices in ([T^ or p^ form a perfect square. For finite 
shear rates wavevector advection leads to a de-phasing of the vertices, resulting in a 
reduction of memory effects. This is the mechanism by which the competition between 
slow structural relaxation and flow enter the theory. Equations P^ - (|14p form a closed 
theory for the calculation of G{t). The corresponding shear stress follows from pT|) . 
For the purposes of this study we consider a one-component system of hard- 
spheres of diameter d, and all numerical calculations will be performed using the 
Percus-Yevick approximation for Sq] discretization of the integrals is performed like in 



Refs. [22|, |49|. This minimal model provides both a reasonable flrst approximation to 
the PMMA colloids used in our experiments and captures the fundamental excluded 
volume effects responsible for the phenomenology observed in our binary mixture 
simulations. The relevant parameters in the hard-sphere system are the bare Peclet 
number Pcq ~ jd'^/Do, measuring the importance of shear relative to Brownian 
motion, and the packing fraction cf) which determines the timescale of structural 
relaxation in equilibrium. The effect of shear on the structural relaxation is measured 
by the dressed Peclet or Weissenberg number Pe, which we define by Pe — jd'^ /{2D), 
using D, the long time self diffusion coefficient in equilibrivim. It compares the shear 



From Equilibrium to Steady State 



10 



10 



10 



,10- 



10 



10' 



0.0 



. 1 1 . 1 1 1 1 . 1 1 . 


. 1 1 . 1 1 1 1 1 1 1 1 
B particles , 


; \x^ 


T=0.14 - 


: Nv ♦ 

- o — o equilibrium 


\"*4 = 


E — --Y=1.2x10"^ 


\ ' - 


: •— -7=3.0x10"^ 


\ ^ 


, ' --'7=6.0x10"'* 


\| : 


- 1 1 1 1 1 1 1 1 1 1 1 1 


' ^ 



2.0 



4.0 6.0 

1/7 



8.0 



10.0 



Figure 3. Self-diffusion constants for the B particles in equilibrium and at the 
indicated shear rates, as obtained from the simulation. For the sheared systems 
the diffusion coefficients are calculated from the mean squared displacement in 
X and y direction, i.e. perpendicular to the shear direction. The temperature 
T = 0.14 is marked. 



rate with the time required for a particle to diffuse a distance of its diameter in one 
spatial direction. 



3. Results 

Figure [3] shows the diffusion coefficients D{j,T) obtained from the simulations of the 
binary Yukawa mixture at various temperatures and shear rates. Only the diffusivity 
of the larger particles is shown, as the smaller particles behave qualitatively identically. 
Open circles represent values measured in the quiescent equilibrium state, D{'y — 0,T). 
They show a strong faster-than-Arrhenius slowing down with decreasing temperature, 
signifying that the system approaches a glass transition. On the timescale of the 
simulation, no diffusive regime is reached in the MSD for temperatures lower than 
r = 0.14. 

Applying steady shear to the system, one finds instead that even at the lowest 
temperatures considered here, the system shows a finite diffusivity. For a fixed value 
of 7, the D(7,T) curves approach a plateau as T is lowered. This plateau occurs at 
values increasingly larger than the quiescent diffusion coefficient. Such an effect is a 
manifestation that the system shows shear-thinning behaviour. At high temperatures, 
a finite 7 has only a small effect on the diffusion coefficient, and consequently all curves 
shown in the figure approach each other for small l/T. 

To study the effects of switching on shear flow suddenly, we consider a fixed 
temperature T = 0.14, corresponding to the lowest temperature where we could 
equilibrate the quiescent system. This allows us to start from a well-defined reference 
state at i = 0, and enables direct comparison with the theory. It also precludes aging 
effects arising from the unsheared state. 

The macroscopic stress <7{t) measured from the time of switching on steady shear 
is shown in Fig. [4] for three different shear rates. The kinetic contribution is negligible 
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Figure 4. Stress-strain relation from the simulation at the temperature T = 0.14 
for three different shear rates, as indicated. 



compared to the potential one, the one studied in the theory. AU three curves are 
quaUtatively similar, and show three regimes. At early times, corresponding to small 
strains 'y = jt < 0.1, (7(7) increases almost linearly with 7. This is the regime of 
solid-like elastic response, and from the prefactor of this almost linear increase one 
can infer the elastic constant. It depends weakly on 7. At large strains, 7^1, one 
finds on the other hand that a{t) approaches an asymptotic value attained in steady 
state, the dynamic shear stress cr(oo) resulting from keeping the system flowing at 
a fixed rate. This corresponds to the regime of plastic deformation. The long-time 
plateau increases with increasing 7, indicating that faster flow induces higher internal 
stresses. The dependence of cr(oo) on shear rate is called 'flow curve'. 

At strains in the intermediate range, 0.1 < 7 < 1, a marked overshoot is seen in 
17(7): the almost linear increase for small strains continues beyond the steady-state 
value of (7, so that a maximum cTraax > cr(oo) results at roughly 10% strain in units of 
the particle diameters, slightly increasing with increasing shear rate. Only after this 
maximum does 17(7) decay towards its steady-state value. Similar features have been 
observed in experiment [ll|, [ij, |l3| and simulation [y, ll6|, |l8|, ll9[ . 



Since the simulation introduces shear flow only through the boundary conditions, 
the sudden switching on of shear does not induce a linear flow profile instantaneously 
throughout the simulation box. Rather, there exists some time scale over which the 
steady-state velocity profile builds up. Especially since we intend to compare the 
simulation results with MCT, where a fixed velocity profile is presumed to exist always, 
let us point out that none of the prominent features shown in Fig. |4] depend on this 
gradual self-adjustment of the profile. To this end, we determined the instantaneous 
velocity profiles, i.e., fs,x(y), for the 7 = 3 x 10~^ simulation at various times 
following shear startup. They have been obtained by averaging over a negligibly 
small time window and over the 250 independent simulation runs. As Fig. [5] shows, 
the profile for t = 0.6 indeed still shows a pronounced 'S'-shaped form, as the Lees- 
Edwards boundary conditions essentially introduce small sheared boundary zones to 
the unperturbed bulk. The propagation of these shear zones inwards is, however, fast. 
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Figure 5. Evolution of the shear profile from equilibrium to steady state, 
as obtained from the simulation of the binary Yukawa mixture. Between the 
equilibrium and the steady state profiles are shown at the times t = 0.6, t = 2.3 
and t = 8.1. The applied shear rate is 7 = 3 X 10~^. The temperature is T = 0.14. 
The inset shows the stress as a function of time for 7 = 3 X 10^'^. The symbols 
in the inset indicate the times along the stress curve at which the shear profiles 
are shown. 



Already around t = 8.1, corresponding to a strain of 7 w 0.02, one essentially recovers 
a linear velocity profile across the whole system within the error bars. This observation 
also holds for larger times (not shown). For comparison, also the profile obtained in 
the steady-state part of the simulation (averaged over a longer time interval and hence 
displaying smaller fluctuations) is shown in the figure. We do not find the formation of 
shear bands or similar localized features in the flow profile, even though the simulated 
equations of motion do not bias towards a linear profile by construction. 

Comparing the strain at which the steady-state flow profile is first attained with 
the typical strains identified from the overshoot features in Fig. [U we conclude that 
the latter are not connected to the peculiarities of whether the system is sheared only 
at the boundaries or forced to follow a linear profile mediated through the solvent 
immediately. This makes comparisons to MCT viable. 

More information on the transient dynamics of the system can be obtained from 
two-time correlation functions, formed between two times t^ > and t > t^^. As 
a simple example, we show in Fig. [S] the mean-squared displacements (MSD) of B 
particles, for three different waiting times, t^ = 0, t^ = 28.1, and t^ sufficiently large 
so that the result as a function of i — t^ becomes tw-independent. For comparison, 
also the equilibrium MSD is shown. It agrees with the MSDs under shear for short 
times. Note that the solid lines correspond to the MSD evaluated in the gradient 
direction, while the dotted equilibrium and steady-state lines represent averages over 
both directions perpendicular to shear. Comparing steady-state and equilibrium, 
we recover for long times the shear-thinning effect discussed in connection with the 
diffusion coefficients, Fig.[3l The equilibrium curve shows a pronounced plateau where 
(r^(f,tw)) ~ 0.013 = Tpjg^^, indicating that particles are at intermediate times caged, 
with a localization length around r' = r^^^^/G « 0.02. This is the usual glassy dynamics 
described by MCT in quiescent systems. The steady-state curve does not show such 
a pronounced plateau, due to the speed-up of the final relaxation. 

As Fig. [6] shows, the sudden start-up of shear flow has a drastic effect on 
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Figure 6. a) Log-log plot of the mean squared displacements of B particles at 
the temperature T = 0.14 as a function of strain 7 = ■yt. Dashed and dotted 
lines correspond to equilibrium and steady state rescaled with 7 = 3 X 10"'^. 
The solid and dashed lines show the transition from equilibrium to steady state 
at this 7 for the vorticity direction {z, 'vd') and the gradient direction {y, 'gd'), 
respectively. Data for the waiting times 4w = and 28.1 are shown (from right 
to left). The dotted lines correspond to the mean squared displacements in 
the equilibrium (EQ) and in the steady state (SS), as indicated. In the inset, 
(r^(i, tw))/(t — iw) is plotted as a function of strain, b) Effective exponent 
/i(t) = d[log(r (t,tyj))]/d[log(t — tw)] for some of the curves shown in a). 
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the MSD measured immediately at the start-up time. While the curve follows the 
quiescent MSD for strains 7 < 0.1, for larger time it suddenly starts to increase much 
more rapidly, and coincides with the steady-state MSD already for 7 w 1, although 
D{'j) ^ D{j = 0). This effect is seen in both directions perpendicular to the flow, 
albeit somewhat less pronounced in the vorticity direction, as demonstrated by the 
solid lines in Fig. \6\ There is thus an intermediate time window where the MSD shows 
super-diffusive behaviour, i.e., grows faster than t. This is even more clearly seen 
when plotting {r'^{t,tw))/{t — iw), as shown in the inset of Fig.[6l At long times, this 
curve would monotonically fall to 2D in the quiescent system. A dip followed by an 
increase at intermediate times, as seen in our data, is the signature of superdiffusion. 

The effect can be quantified more precisely through a logarithmic derivative, 
/z(i) = (i[log(r^(i,tw))]/d[log(i — iw)], which approaches ^ = 1 for ordinary diffusion 
and /i = 2 for ballistic motion. The simulation result for /i(i) and iw = is shown 
both for the vorticity direction and the gradient direction in the lower panel of Fig. [HI 
At short times, ^{t -^ 0) ^ 2 reveals the ballistic motion underlying the molecular 
dynamics simulation. At t w 1, fi(t) « is indicative of the equilibrium cage effect. 
In the iw = curve, it is truncated at jt w 0.1, and quickly fi{t) raises to about 
2 in the gradient direction, and 1.8 in the vorticity direction, before finally settling 
to /U,(i -^ 00) -^ 1 as expected. The simulation thus reveals an almost 'ballistic' 
regime in the MSD at intermediate times. Remarkably, this result remains unchanged 
if one performs the simulation with a much larger damping coefficient in the DPD 
thermostat, ^ = 1200 [33|. This highlights that the regime of superdiffusion is likely 
a feature of sheared glassy dynamics and not of the short-time microscopic motion. 
It is worth noting that /x < 1 holds for all times beyond the microscopic transient in 
the steady-state curve. This is what one expects from the regime of true 'structural 
relaxation', and what holds precisely in the Brownian equilibrium system, where all 
correlation functions are purely relaxing functions. 

Figure [7] displays the non-Gaussian parameter a2(i), as defined by Eq. ^, for 
equihbrium, steady state and waiting time iw = at T = 0.14 and 7 = 3 x 10~^. Also 
shown is the equilibrium a2{t) at T = 0.18 where D^j = 0) is approximately equal to 
£'(7 = 3 X 10^'^) at T = 0.14. In the steady state the amplitude of maximum in a2{t) is 
about an order of magnitude smaller than that for the equilibrium at T = 0.14 and also 
slightly smaller than that of the equilibrium curve at T = 0.18. Thus the application 
of shear seems to reduce non-Gaussian effects at intermediate times. Similar to the 
behaviour of the MSD, the non-Gaussian parameter for waiting time iw = follows 
first the equilibrium curve. Then, at time t > 30 the decay of a2 can be fitted to a 
power law, ^2 oc t~"-^'^. 

Let us now compare in detail the findings from the simulation to the colloidal 
experiment. Figure [8] shows the mean-squared displacements obtained from confocal 
microscopy and particle tracking, taken in the vorticity direction. Here only a single 
shear rate is shown. Note in comparing with Fig. [HI that the relevant quantities 
characterizing the magnitude of shear differ between the two setups: while the 
experiment corresponds to a Weissenberg number Pe « 7 (Pcq « 1), it is closer 
to Pe « 288 in the simulated system. Apart from obvious differences due to this 
circumstance, the results are qualitatively similar. As the inset of Fig. [8] demonstrates 
for the experimental {r'^{t, tv,)) / {t — t^) , again a nonmonotonic dip is found at roughly 
7 w 0.1, identifying a corresponding super-diffusive regime in the MSD. Like in the 
simulation, this effect is most pronounced for t^ ~ and then continuously weakens 
as the steady state is approached. 
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Figure 7. Non-Gaussian parameter 0:2 (f) for equilibrium, steady state and 
waiting time iw = at the temperature T = 0.14. For the sheared case, the 
shear rate is 7 = 3 X lO" ^. For comparison also the equilibrium a2{t) at T = 0.18 
is shown. _D(7 = 0) at this temperature is approximately equal to D(7 = 3 X 10"^) 
at T = 0.14. The bold solid line is a fit to the tw = curve with the indicated 
power law 02 tx i^^.SS -pj^g inset shows experimental results for a2{t). Here, the 
shear rate is 7 = 4.5 X lO".^ s~^. 
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Figure 8. Log-log plot of the experimental mean squared displacement in 
vorticity direction as a function of strain ■yt for different waiting times, as 
indicated. The inset shows {r^{t, t^)) /{t — iw) as a function of 7-4. 
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Figure 9. Steady state shear profile, as obtained from the experiment (dots). 
Dashed lino: linear fit providing 7 = 0.045 s~^. 



A quantitative difference is found when comparing the extension and strength 
of the super-diffusive window in the MSD. In the simulation, the crossover in the 
iw = curve from the equihbrium to the steady-state hmiting cases is much more 
rapid. Likewise, the experiments show the first increase of the transient MSD beyond 
the equihbrium curve for shghtly smaller strains, 7 « 0.05, than the 7 ~ 0.1 found in 
the simulation. 

Also the non-Gaussian parameters have been determined in the experiment. They 
are shown in the inset of Fig. [T) One can infer that also this quantity behaves 
qualitatively as in the simulation. 

Again it is important to stress that the effects discussed here are not consequences 
of shear localization or, more generally, developments of non-linear flow profiles. This 
has also been checked for the experiment. The velocity profile in steady state is almost 
perfectly linear, as Fig. [9] demonstrates. This confirms that also in the experimental 
setup, a homogeneous constant-shear steady state is approached, without obvious 
shear-localization features or wall effects. 

Having established the close connection between simulation and experiment, let 
us now try to rationalize these findings with MCT. The quantity a{t) is easily accessed 
in the mode-coupling framework and thus forms a convenient starting point for the 
analysis. In order to provide a meaningful comparison with simulation results, we 
have to first map the relevant quantities between the theoretical hard-sphere model 
and the binary Yukawa mixture employed in the simulations. Since we are dealing 
with slow structural dynamics, it is natural to fix the final relaxation time r first 
for equilibrium and then to consider the effect of shear identifying the Peclet and 
Weissenberg numbers. The only free parameter in the quiescent hard-sphere model is 
the packing fraction, adjusted in the following to Aip = (j> — cf>c = —1.16 x 10~^ from 
the requirement that MCT describes the equilibrium MSD found in the simulation for 
long times (see below). Next, the shear rate appropriate for use in the theory may be 
determined. We find Pe = 183 (corresponding to 7 = 5.5 x 10~^Do/(P) to lead to 
the same D{'y)/D{'j = 0) as it arises from the simulation at Pe w 288. 

In Fig. [TOl we show the resulting MCT-calculated stress a{t) as a function 
of accumulated strain jt for the three shear rates considered in the simulation. 
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Figure 10. Comparison of the stress from simulation (dots) and theory 
(solid lines) for the shear rates considered in simulation, 7 = 3 X 10~^ (green), 
7 = 6 X 10~* (red) and 7 = 3 x 10~* (black). The relation between simulation 
energy units (-Enid) ^-i^d theory is -Bnid = 8.8kgT. The theoretical shear rates 
(in units of d? /Dq) obtained from fitting the MSD at the highest shear rates of 
the simulation (corresponding to a theoretical value Pc=183) are 7 = 5.5 X 10"'^, 
7 = 1.1 X 10~^ and 7 = 5.5 X 10"'*, respectively. 



Theoretical results are shown in hard-sphere units [ksT /d?) whereas simulation results 
are expressed in units of i?md = t/d^. The two axes have been adjusted in order to 
match the stress attained in steady state, fixing a relation of i?md ~ S.SfcsT between 
the kinetic energies of the different systems for the largest shear rate considered in 
the simulation. This fitted value of T is surprisingly close to the actual value in 
the simulation, indicating that the theory captures the stress magnitude well. As 
Fig. [To] shows, this also provides a reasonable description of the initial small-strain 
stress, reflecting the elastic constant of the system. More importantly, the theory 
qualitatively reproduces the stress overshoot found in the simulation. That the peak 
occurs at strains of roughly 7 « 0.1, in accord with simulation, is a consequence of our 
choice for the parameter 7c in the isotropic approximation to the advected wavevector. 
The magnitude of the overshoot is significantly underestimated, which might be due 
to the MCT approximation per se, but may also be a result of the additional isotropic 
approximation underlying our calculations. In the decay of a{t) after its maximum, 
the simulation data show a rather slow decay for log]^Q7 > 0.75, which is also not 
captured in MCT. 

The theory quite naturally relates this stress overshoot to a peculiar feature of the 
dynamic shear modulus G{t). In Fig.llH this modulus is shown both in the quiescent 
case and for the shear rates used in Fig. [101 The equilibrium function shows the 
two-step relaxation pattern familiar from other correlation functions close to the glass 
transition. The non-linear shear thinning is clearly seen as a speeding up of the final 
relaxation time in G{t) from the equilibrium case (r « 10"^) to that of steady-shear, 
for which the characteristic time scale is 7^^. A plateau at times around t « 1 is still 
seen as a remnant of the cage effect. Note that only much closer to the glass transition 
would a well-established constant plateau develop, defining the solid shear modulus 
Goo. For the parameters chosen here, the modulus G{t) still shows a finite slope in 
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Figure 11. Generalized shear modulus from mode-coupling theory, calculated 
using ij) ~ 4'c = —1.16 X 10"'^, for 7 = (rightmost curve) and 7 = 5.5 X 10~'^, 
1.1 X 10"'^, and 5.5 X 10~* (left to right), corresponding to curves in figure [TUl 
The stress overshoot in Fig. 1101 is related to the negative undershoot in G(i) at 
long times. The plateau value which would develop closer to the glass transition 
is indicated by the broken line. The inset shows the moduli as a function of strain, 
where the vertical bars indicate the strains considered in Fig. 1121 




Figure 12. The product S'^S',,.-. which enters the stress vertex and which 

leads to negative values of the modulus G(i) at long times. The uppermost panel 
and inset show S',S',,.^ and S^, respectively, for a packing fraction (f> = 0.4. 

Curves are for strain values -yt = 0, 0.1, 0.2, 0.3, 0.4 and 0.5 (where the arrow 
indicates the direction of increase). The lower panel shows the same quantities 
for a packing fraction of </< = 0.5148, corresponding to the packing fraction used 
in our comparison with simulation. The five values of the strain considered are 
marked in Fig. 1111 



the plateau regime. Considering that a{t) is the time-integral over G{t) leads to the 
conclusion that the 'elastic' small-strain regime visible in Figs. |4] and [10] is not truly 
linear in 7. Indeed, the simulation data could be fitted with a power law cr oc 7^ with 
some cfFectivc exponent a; < 1 in this time window. 
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The transient shear modulus in addition shows a charaeteristic negative dip for 
strains around jt w 0.1, which for the highest shear rate corresponds to IDq/cP w 20. 
From the integral for the stress, Eq. (fTTI) . it is apparent that an overshoot in (j{t) 
is the result of such a negative portion in G{t) for long times. Recall that G{t) is a 
stress-stress autocorrelation function, so that this negative dip may be interpreted 
as a 'backscattering' of stresses. It arises from the final stage of the relaxation, 
commonly referred to as the a process, which is attributed to the breaking up of 
nearest-neighbour cages. As a result MCT does not correlate the stress overshoot to 
peculiar inhomogeneities, such as localized shear events. 

Within our MCT-based approach negative values of the modulus occur as a result 
of the factor S'l^SL^-. appearing in Ea. (fT2)) . The derivative S'^. oscillates about zero as 
a function of k and attains maximal positive and negative values in the vicinity of 
the first peak in Sk- As the strain is increased the resulting de-phasing of the two 
factors S'l. and S'^s leads to negative regions in the product S'^.SL^y This behaviour 
is shown in more detail in Fig[T2] for two values of the packing fraction (j), where the 
structure factors are taken from the analytic Percus-Yevick theory. Comparing the 
curves for the two densities it is apparent that the first peak in Sk is considerably 
narrower and steeper for (p = 0.5148 than for (f) — 0.4, leading to larger values of 
S'^ in this region. In each case the product S'^.S'Lj) becomes increasingly negative as 
the strain is increased from ^t = 0.1 to ^t = 0.2. At longer times, the decay of the 
correlators in Eq. p2)) suppresses the negative vertex. The inset to Fig. [11] shows the 
modulus G(i) as a function of strain. For strains in the range 0.1 to 0.25 the modulus 
attains its most negative values. Comparison of the two panels in Fig. [12] shows that 
for = 0.5148 the negative region in S'/^S'^, is much more developed at jt = 0.2, 
relative to = 0.4. This reflects the differences in the first peak of the structure 
factor at the two densities and suggests that systems with a steeper and narrower first 
peak in Sk will possess a more pronounced overshoot within the present MCT-based 
approach. Structure factors with a steeper and narrower first peak are known to occur 
in two-dimensional systems [50| . 

We have now established two features of the transient dynamics immediately 
following the instantaneous switch-on of shear, both of which arise at strains of about 
10% of the particle diameter. The first of these, a{t) (equivalently G(i)), concerns the 
collective dynamics, while the other, {r'^[t,t^)), concerns the single-particle dynamics 
of a tracer particle. It is intriguing to see in which way these conceptually different 
pieces of information might be related and to connect the macroscopic rheological 
information contained in a{t) with the microscopic particle motion. We consider in 
more detail the transient MSD, 6r'^{t) = {r'^{t,t^ ~ 0)). An idea that has proven 
particularly successful in the quiescent system is that of a generalized Stokes- Einstein 
(GSE) relation, proposed by Mason and Weitz as a correspondence between macro- 
and micro-rheology [51j . Within this approach the diffusion coefficient D is related to 
the viscosity r] by 

ksT ksT 

indarj Sirda J^ ds G{s) 

In the context of q uiescent MCT, this relation exemplifies the coupling of the a- 
relaxation times |52l | , an d even quantitatively works surprisingly well for the quiescent 
hard-sphere system [49[ . Quantitatively accurate results can be obtained if one allows 
for a small correction to the geometrical prefactor Sird appearing in the original Stokes- 
Einstein relation, a ^ I. 
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Figure 13. Comparison of the mean-square displacements from simulation 
(vorticity direction, dots) and theory (solid lines), for equilibrium (EQ) and for 
the highest shear rate considered in the simulations (waiting time tw = 0). A 
superdiffusive regime is apparent in the theoretical MSD for t^ = 0, establishing 
the connection to the stress overshoot in Fig. 10. Inset (a) shows the effective 
exponent from the theory for tw = obtained from the logarithmic derivative. The 
peak is at a value of ~ 1.2. Inset (b) shows the MSD divided by t. Dash-dotted 
lines correspond to the MSD calculated from the simulated cr{t) at 7 = 3 X 10"'^ 
and the GSE relation (see text). The dashed line indicates the linear diffusive 
behavior in steady state. 



In addition to the correlator $fc(t), MCT also yields equations of motion for the 
tagged-particle correlation functions, which yield an equation for the MSD in the q — > 
limit. Taking the form of the equation of motion for Sr'^{t) from the quiescent theory 
(ignoring any generalizations that come about due to the anisotropic flow geometry), 
we can write 

^r^ {t) + ^ f dt' 7«(^) (t - t')5r'' {t') 



keT 



6Dot, 



(16) 



where m^^'>{t) is the tagged-particle memory kernel evaluated as q ^ 0. An analysis 
of the above equation at large times immediately gives an expression for the diffusion 
coefficient in terms of the memory kernel, 

D= ^ « I (17) 

Comparison of ([T7]) with psp clearly establishes a relationship between the time 
integrals of to'*' (t) and G{t). However, the connection is deeper than simple equality of 
the integrated quantities. Studies performed using the quiescent mode coupling theory 
have established that the MCT approximation to G{t) is in excellent agreement (up 
to a constant prefactor) with the MCT approximation to m*^*' (i) for the hard sphere 
system [49[ . In this work we make the assumption that this correspondence holds also 
in the non-linear regime for which G{t) = G(t,'^). Given this assumption we obtain 
a direct correspondence between the MSD and the derivative of the shear stress after 
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startup, {d/dt)a{t) ex G{t), 

We now consider application of the GSE approximation (fT8|) to the calculation 
of the mean-squared displacement using ([T6|) . In Fig. [13] we compare the tw = 
curves obtained from simulation and theory for the largest shear rate considered in 
the simulation. Also shown is the equilibrium curve that has been used to adjust 
A(j) in fitting the equilibrium simulation result. At short times differences arise due 
to the different microscopic dynamics, as expected, but for the long-times relevant 
to our discussion, the fit is convincing. The MCT-GSE theory displays a super- 
diffusive regime at intermediate times, qualitatively reproducing the phenomenology 
of the simulation. In order to quantify this superdiffusivity, we show in inset (a) the 
logarithmic derivative of the MSD, yielding an effective exponent for 5r'^{t). While 
in the simulation this quantity grows up to 2 (reflecting ballistic motion) in the 
intermediate regime, our MCT calculation produces effective exponents up to roughly 
1.2. This situation is much closer to the results of our colloidal experiments. 

Note that from ([T7]), a super-diffusive regime in the MSD can only arise if 
m(*^(t) takes on negative values. By virtue of pS]) . the presence of a super-diffusive 
regime is therefore intimately connected to the stress overshoot and the corresponding 
overrelaxation feature found in G(t). It is therefore tempting to assume that the more 
pronounced super-diffusive regime found in the simulations relative to the theory 
is a consequence of the correspondingly stronger stress overshoot (see Fig. [TU|) . In 
order to investigate this connection further we have directly calculated the mean- 
squared displacement from the simulation curve for cr(t), using (jlSp to obtain the 
memory function m^'''{t). In this calculation, the Newtonian-dynamics equivalent of 
Eq. (|16p was used with a prefactor a = 0.833 (accounting for the numerical uncertainty 
that arises from the neccessary smoothening of the derivative of the simulated (j{t)). 
The resulting MSD is shown in Fig. [T3] (dash-dotted lines). For times longer than 
tDo/d^ sa 100 the MSD obtained in this fashion closely follows the directly simulated 
quantity. However, for shorter times the agreement is less satisfactory and significant 
deviations are visible. The logarithmic derivative shown in inset of Fig. 113b confirms 
this. As a result of these findings we conclude that the approximation psp remains 
qualitatively reliable in the non linear regime but that it fails quantitatively in the 
super-diffusive time window where the stress overshoot is strongly pronounced. 

Finally, we consider application of the GSE-MCT relation to fitting the 
experimental data for the MSD, the results of which can be seen in Fig. [TH We 
follow the same procedure as when fitting the simulation data and consider first the 
equilibrium MSD. By adjusting Acj) wc fix the ratio of final to microscopic relaxation 
time and find a value A<j> = —1.1 x 10^^ reproduces the experimental results to a 
satisfactory level, using also a value a = 0.95 little different from unity. We next 
increase the shear rate until the long time diffusion matches that of the simulation. 
This is achieved for Pe = 16 (corresponding to 7 = 5 x 10^^ Dq / d^) . We find a super- 
diffusive regime in qualitative agreement with the experiment at Pe = 7, but which is 
somewhat underestimated in magnitude. This lends support to our argument that the 
approximation (|18p is less reliable in the super-diffusive regime than in equilibrium. 
In the inset wc show the MSD divided by t which confirms the relatively weaker 
superdiffusion. 
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Figure 14. Theoretical fits to tlie experimental mean-square displacement 
(dotted) for equilibrium and iw = (solid lines) using a = 0.95 (see text). The 
inset shows the MSD divided by t and highlights the super-diffusive regime. 



4. Conclusions 



We have studied the transient dynamics of a glassforming system that is subjected 
to the sudden commence of steady shear flow. Results from computer simulation, 
confocal-microscopy experiment, and mode-coupling theory give a consistent picture 
that reveals peculiar dynamical features around the time where the system reaches 
10% strain for the first time. In particular, the mean-squared displacements in both 
the vorticity and the gradient direction show super-diffusive behaviour at intermediate 
times, interpolating between the equilibrium curve that is followed for 7 ;$ 0.1 and 
the faster steady-state curve that is followed for 7 > 1. The super-diffusive behaviour 
is more pronounced in the (ballistic) simulation, giving rise to effective exponents of 
2, and less pronounced in the colloidal experiment. 

Again at strains of 7 ~ 0.1, a well-recognized feature in the collective dynamics 
is found, visible in the building up of the shear stress as a 'stress overshoot' (a 
maximum in a{t) followed by a decrease towards a smaller steady-state value at later 
times). While we cannot propose a simple picture of microscopic motion behind 
this stress overshoot, we are able within MCT to directly relate it to a 'stress 
overrelaxation', visible in the corresponding stress-stress autocorrelation function in 
the late a-relaxation regime. It can thus be attributed to the peculiar way in which 
cages break up, and MCT predicts its specific shape and magnitude to be system- 
dependent. We have demonstrated that MCT in this way qualitatively describes the 
stress overshoot seen in the simulation, albeit quantitatively underestimated. 

A generalized Stokes-Einstein relation can be used to relate the above two features 
of the transient collective and single-particle dynamics. We have used MCT in 
combination with this close link between micro- and macro-rhcology in order to 
demonstrate that the theory-predicted stress overrelaxation quite naturally leads to a 
super-diffusive regime in the mean-squared displacement. Again, its magnitude is less 
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pronounced than in the simulation, but compares reasonably well with experiment, so 
that the remaining differences can be attributed to the microscopic difference in the 
systems, the isotropic nature of our approximations, and violations of the generalized 
Stokes-Einstein relation. 
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